Astronomy & Astrophysics manuscript no. Suad-astroph 


©ESO 2011 


December 5, 201 1 





GS1 00-02-41 : a new large H I shell in the outer part of the Galaxy. 

L. A. Suad 1 , S. Cichowolski 2 , E. M. Arnal 1 - 3 , and J. C. Testori 1 

1 Institute Argentines de Radioastronomfa (IAR)(CCT - La PLata), CC 5, 1894, Villa Elisa, Argentina. 

2 Instituto de Astronorma y Fi'sica del Espacio (IAFE), CC 67, Sue. 28, 1428 Buenos Aires, Argentina. 

3 Facultad de Ciencias Astronomicas y Geofisicas, Universidad Nacional de La Plata, Argentina. 

received /accepted 

ABSTRACT 

Context. Massive stars have a profound effect on the surrounding interstellar medium. They ionize and heat the neutral gas, and due 
to their strong winds, they swept the gas up forming large HI shells. In this way, they generate a dense shell where the physical 
conditions for the formation of new stars are given. 

Aims. The aim of this study is to analyze the origin and evolution of the large HI shell GS 100-02-41 and its role in triggering star 
forming processes. 

Methods. To characterize the shell and its environs, we carry out a multi-wavelength study. We analyze the HI 21 cm line, the radio 
continuum, and infrared emission distributions. 

Results. The analysis of the HI data shows an expanding shell structure centred at (/, b) = (10O;6, -2°04) in the velocity range from 
-29 to -51.7 kms~'. Taking into account non circular motions, we infer for GS 100-02^1 a kinematical distance of 2.8 ± 0.6 kpc. 
Several massive stars belonging to Cep OB 1 are located in projection within the large H I shell boundaries. The analysis of the radio 
continuum and infrared data reveal that there is no continuum counterpart of the H I shell. On the other hand, three slightly extended 
radio continuum sources are observed in projection onto the dense HI shell. From their flux density determinations we infer that 
they are thermal in nature. An analysis of the H I emission distribution in the environs of these sources shows, for each of them, a 
region of low emissivity having a good morphological correlation with the ionized gas in a velocity range similar to the one where 
GS 100-02-41 is detected. 

Conclusions. Based on an energetic analysis, we conclude that the origin of GS 100-02-41 could have been mainly due to the action 
of the Cep OBI massive stars located inside the HI shell. The obtained age difference between the HI shell and the HII regions, 
together with their relative location, led us to conclude that the ionizing stars could have been created as a consequence of the shell 
evolution. 

Key words. ISM: structure - ISM: kinematics and dynamics - HII regions - Stars: formation 



1. Introduction 

The presence in the interstellar medium (ISM) of the Milky Way, 
when viewed in the A ~ 21 cm line emission of the neutral hy- 
drogen (H I), of giant structures having linear dimensions of a 
few hundred parsecs in diameter is a widely kno wn and well ob- 
served phenomena first noticed by iHeilesf dl979b . These struc- 
tures are usually detected as huge shells or arc-like features of 
enhanced H I emission surrounding regions of low H I emissiv- 
ity, receiving the generic name of H I supershells. These features 
may even be the dominant structure in the interstellar medium, 
taking up a large fraction of the volume of the galactic disk. The 
HI structures could be also observed at infrare d wavele ngths. 
Based on the 60 and 100 pm IRAS databases, iKonvves et alj 
(2007) have performed an all-sky survey of loop- and arc-like 
structures. 

Similar st ructures have also been observed in nearb y spi- 
ral galaxies riStanimirovicll2007t IChakraborti & Ravi 1201 lT) . In 
the Mi l ky W a y, the se structures were initially catalogued by 
iHeilesi (119791 Il984h . Though a large number of HI features 
likely to be classified as either large HI shells or HI su- 
pershells have been catalo gued in the outer (90° < I < 
270°) part of the Galaxy (lEhlerova & Palousl 12005b . only a 
small number of them have been studied in some detail. The 
later dJung et all 1 1 9961: IStil & Irwir] 1200 ll lUvaniker & Kothesl 
l2002HMcClure-Griffifhs et al.l2002l:ICazzolato & Pineauld2 003: 



Arnal & Corti 2007; Cichowolski & Pineault 201 1) have galac- 
tocentric distances ranging from 9.7 to 16.6 kpc, diameters from 
120 to ~ 840 pc, expansion velocities between ~10 and ~ 20 
kms~', and kinetic energies from ~ 1 x 10 50 up to ~ 6 x 10 51 
erg. A mong them only two (GS 305 +01-24 and the feature stud- 
ied by ICazzolato & Pineaultl d2003b ) have an OB-association as 
their likely powering source, and other three (GSH 9 1 .5+2-1 14, 
GS 234-02 and GS 263-02+45) show evidence of having in- 
duced the formation of new generation of stars. 

The general consensus is that those structures whose ki- 
netic energy is of the order of, or less than, a few times 10 51 
erg, very likely may have been created by the joint action of 
stellar winds and supernova explosions. A large number of ex- 
amples are reported in the literature ("e.g. |Uvaniker & Kothesl 
120021: lArnal & Cortil 120071: ICichowolski & Pineaultl 1201 lb . On 
the other hand, for expanding H I structures having kinetic en- 
ergies in excess ~ 10 52 erg, termed supershells, the above 
mechanism may not be adequate to create them because one 
would need a stellar grouping (either an open cluster or an OB- 
association) with many more stars than the average found in the 
Milky Way. In these cases alt ernative mechanisms like the in- 
falling of high velocity clouds dTenorio-Taglelll98lb or gamma- 
ray bursts (Perna & Raymond 2000) may be at work. 

Along its expansion, these structures (either a shell or a su- 
pershell) may became gravitationally unstable, forming clouds 
that later on may lead to the formation of stars along the periph- 
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ery of these H I structures, or else the expanding structure may 
hit and compress pre-existing ISM molecular clouds from one 
side. During this process a high density perturbation may move 
into the molecular clouds, which may eventually collapse into 
denser cores in which star formation may occur. A thorough re- 
view of ob servations and theo ry related to trigger star formation 
is given by Elmegreenl (Il998l) . 

A new large scale study aim at detecting in the outer part of 
the galaxy structures likely to be either large H I shells or super- 
shells is being carried out by one of the authors (L. A. Suad) as 
part of her PhD Thesis. 

In this paper we analyze a new large H I shell observed at (I, 
b) ~ (100°, -2°), with the purpose of elucidating both its origin 
and its interaction with the surrounding ISM. We also look for 
signs of recent star formation activity likely to be related to this 
shell. 



Table 1. Observational parameters. 



2. Observations 

Low resolution HI data were retrieved 
Argentine-Bonn (LAB) survey 



from the L eiden- 
dKalberla etail 120051) . This 
database is well suited for a study of large scale structures due to 
its angular resolution. The entire database has been corrected for 
stray radiation. High angular resolution H I data, covering most 
of the structure under study, w ere obtained from the Canadian 
Galactic Plane Survey (CGPS. iTavlor et al. 1 120031) . Besides the 
high resolution HI data, the CGPS also pro vides high reso- 
lution continuum data at 408 and 1420 MHz ( Landecker et al J 
2000) . Continuu m data at 2695 MHz (iReich et alJll984 11990: 
Fiirst et al 1990) were also used in this study. 

Infrared data f rom the Midcourse Space Experiment (MSX) 
dPrice et al] 12001 1) were obtained from the Infrared Science 
Archive Q Infrared image s from the Improved Reprocessing of 
the IRAS Survey (IRIS) (iMiville-Deschenes & Lagachd kOOS) 
were retrieved from the Sky View homepageB- 

Carbon monoxide data ( 12 CO (1— »0)) for a small area cen- 
tered at (/, b) ~ (100°7, -0°5) were kindly made available to us 
by Dr. Christ Brunt, while for / > 102°08 12 CO data were re- 
trieved from the CGPS database. These data were obtained using 
the 14-mdish of the Five College Radio Astronomy Observatory 
(FCRAO). 

TableQ]summarizes the most relevant observational parame- 
ters. 

The use of this multi-wavelength approach allows us to 
probe the different components (neutral and ionized gas and 
dust) of the ISM. 



3. GS 100-02-41. 

3.1. Neutral hydrogen data. 

We consider that a given HI structure may be classified as a 
shell if the following criteria are fullfilled: 

1 . It must have a well defined lower brightness temperature sur- 
rounded (partially or completely) by regions of higher tem- 
perature. 

2. The HI minimum must be observable in at least 5 consecu- 
tive velocity channels. 



1 The NASA/IPAC Infrared Science Archive is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, un- 
der contract with the Nasa Aeronautics and Space Administration 
(http://irsa.ipac.caltech.edu I 

2 http://skyview.gsfc.nasa.gov/ 



LAB HI data 




Angular resolution 


34' 


Velocity resolution 


1.3 kms" 1 


Velocity coverage 


-450 to 400 km s' 


CGPS HI data 




Angular resolution 


L2x 1' 


Velocity resolution 


1.3 kms- 1 


Velocity coverage 


-164 to 58 kms" 


Radio continuum 




Angular resolution (408 MHz) 


3:5x416 


Angular resolution (1420 MHz) 


L2X 1' 


Angular resolution (2695 MHz) 


4:3 


CO data 




Angular resolution 


~ 1' 


Velocity resolution 


0.824 kms" 1 


Infrared data 




Angular resolution (MSX) 


18'.'4 


Angular resolution (IRIS) 


3:8-43 


Angular resolution (HIRES) 


0:5 - 2' 



3. It must have a minimum angular size of 2°. This condition is 
set bearing in mind the angular resolution of the LAB survey, 
to assure the structure under study is, angularly speaking, 
fully resolved. 

4. At the kinematic distance of the structure, its linear size must 
exceed 200 pc. 

Relying on the criteria described above, we are constructing 
a large H I shells and supershells catalogue in the outer part of 
the Galaxy (Suad et al. in preparation). As part of this catalogue, 
we have found an HI structure located at (I, b, V) = (100°6, 
-2°3, —41 kms -1 ). Following the standard nomenclature, this 
shell is labelled as GS 100-02-41 . 

To illustrate the main observational findings, in Fig.Q]a mo- 
saic of six HI images covering the velocity range form -47.41 
to -31.96 kms -1 is shown (all velocities in this paper are re- 
ferred to the LSR). Each image is an average of three individual 
line channel maps and covers a velocity interval of about 3.09 
km s _1 . 

The upper left panel shows the H I brightness temperature 
distribution at -47.41 kms -1 . Besides the HI emission arising 
from the galactic plane, quite a few regions of low emissivity are 
clearly identified. These regions appear as a local minimum sur- 
rounded along most of its perimeter by regions of higher bright- 
ness temperature, as for example the structure located at (/, b) 
— (10L0, -2°0). As we move towards more positive velocities, 
this structure becomes a well defined H I minimum within the 
overall galactic H I emission. This minimum achieves its maxi- 
mum angular extent at -41 .2 km s~' and is barely visible at -32 
kms -1 (lower right panel). 

An average of the H I emission distribution in the velocity 
range from -42.26 to -34.02 kms -1 is shown in Fig. [2 where 
a huge HI structure having dimensions of ~ 6° x 3° (Al x Ab) 
is easily recognizable. 

To estimate the main parameters of this large H I shell, we 
have characterized the ellipse that best fit it using a least square 
method. The points used to fit the ellipse correspond to the lo- 
cal maxima around the cavity. To select these points the CGPS 
H I data were used. The parameters derived from this fit are: the 
symmetry centre of the ellipsoidal HI distribution (Iq, bo), the 
length of both the semi-major (a) and semi-minor (b) axes of 
the ellipse, and the inclination angle (6) between the major axis 



2 



L. A. Suad et al.: GS100-02-41: a new large HI shell in the outer part of the Galaxy. 



19 OT -3J « M WO 7D » 




IM" 10^" H* #** 



lDff IMF 10*" UK* W T 9^ M* 



K*.' iD|" TCHf 1 IE" -Ktf BC M" M 1 



f*rn 1 1 — - c — r 




TajL4f$ -3i.fi. 



- -ii 




l«' l«* l<M* w* W *** 



im* so** id*' hi* i<wP w* iff *+*■ 



iwp "uf" ^h" i «j j." '-ac 1 &v <h' 



Fig. 1. LAB HI emission distribution averaged in the velocity range from -47.41 to -31.96 kms~'. Each image is an average of 
three velocity channel maps. Contour levels are at 20, 30, and 40 K. The LSR central velocity of each image is indicated in its top 
left corner. The color bar-code shows the temperature scale in kelvins for all the images. 



and the galactic longitude axis. This angle is positive towards the 
north galactic pole. They are given in Table|2] 

Under the assumption of a symmetric expansion, an ellip- 
soidal HI feature having a central velocity Vo and an expan- 
sion velocity V e should depict, in a position-position diagram, an 
ellipse-like pattern when observed at different radial velocities. 
At Vo the ellipse of H I emission attains its maximum dimen- 
sions, while at extreme velocities (either approaching (V m = 
Vo - V e ) or receding (Vm = Vo + V e )) the hydrogen emis- 
sion should look like an "ovoidal" patch of emission. At inter- 
mediate velocities the dimension of the HI ellipse shrinks as 
Vm (or V m ) is approached. The expansion velocity is estimated 
as half of the total velocity range (V e = 0.5 (| Vm - V m |)) cov- 
ered by the H I emission related to the feature. This method al- 
ways provides a lower limit to V e , because HI emission arising 
from those regions having radial velocities close to either the 
maximum approaching (V m ) or receding (Vm) cap are usually 
difficult to disentangle from the overall galactic H I emission. 
Another way to determine V e without the above drawback is to 
use the velocity-position diagrams. In Fig. [3] a radial velocity 
versus galactic longitude diagram is shown. This image shows a 



different view of GS 100-02^41 . The inner region of the ellip- 
tical HI feature shown in Fig. [2] corresponds to the low emis- 
sivity region seen at (V, I) = (-41 kms -1 , 100?6). The strong 

peaks of HI emission seen at 29 kms~' and 51 kms" 1 

arise from the H I associated with the expanding walls of neutral 
gas related to GS 100-02—41. Figure|4]shows a cross-cut of Fig. 
[3] along / = 100°6. By making a Gaussian fit to the HI emis- 
sion peaks, V m and Vm can be derived, along with the value of 
Vo, that corresponds to the velocity between V m and Vm where 
the minimum value of Th is observed. We obtain V) = -41 ± 2 
kms~'and V e = 11 ± 2 kms -1 . 

It is well known that non-circul ar motions on a larg e scale 
are present in the Perseus spiral arm dBrand & Blitzf l993). mak- 
ing the galactic rotation model not suitable for this region of 
the Galaxy. Based o n the observed velocity field derived by 
iBrand & Blit3 (1 19931) . the systemic velocity of GS 100-02-41 
can be translated to a kinematic distance of 2.8 ± 0.6 kpc. We 
shall adopt this distance for the H I shell. 

Under the assumption that the H I emis sion is optical ly thin 
and following the procedure described bv iPineauIdd 19981) . the 
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Fig. 2. LAB H I emission distribution averaged in the velocity 
range from -42.26 to -34.02 kms -1 . Contour levels are from 
10 to 90 in steps of 10 K. The area below the dashed line cor- 
responds to the region where CGPS data are not available. The 
symbols indicate the location of the star members of CepOBl 
lying inside the ellipse (see Section [3~3l . White star symbols cor- 
respond to O-type main sequence stars. Asterisks and triangles 
indicate evolved stars that were O-type and B-type stars during 
the main sequence phase, respectively. 
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Fig. 3. LAB H I emission distribution averaged in the latitude 
range from -3°0 to -1°3. Contour levels are from 13 to 53 K in 
steps of 5 K. The dashed line indicates the location where the 
cross-cut shown in Fig. [4] was performed. 

total neutral hydrogen mass of a structure located at a distance d 
(kpc) that subtends a solid angle Q. (square arc-min) is given by 

M ra (M ) = 1.3 x 10~ 3 d 2 AVAT B Q. 

where AV is the velocity interval over which the structure is de- 
tected, expressed in kms -1 and ATb (K) is the mean brightness 
temperature defined as AT% =| r s h - 7b g |, where r s h refers to 
the mean brightness temperature of the H I shell, and T\, g corre- 
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Fig. 4. HI cross-cut at 1=100?6 obtained from Fig. [3]. The ap- 
proaching (V m ), systemic (Vo) an d receding (V M ) velocities are 
indicated. 

sponds to the temperature of the contour level defining the outer 
border of the H I shell. The latter represents the temperature of 
the surrounding galactic HI emission gas. For this structure we 
estimated T sh = 33.5 ± 1.6 K, T bg = 21.5 ± 1.5 K, AV = 22 + 2 
kms^and £2 = 4.2 x 10 4 arcmin 2 . Adopting solar abundances, 
the total gaseous mass of GS 100-02^11 is M t (M ) = 1.34 Mm 
(see Table |2]i . 

Assuming that the mass is uniformly distributed within the 
structure's volume, we obtain the gas number density (n s /,) of 
the swept up gas. The ambient density (no) of the medium into 
which the large HI shell is evolving is derived by uniformly 
distributing the shell mass (M t ) over the volume swept up by the 
structure. Assuming this mass to be distributed into an ovoidal 
volume whose semi-major axis is a and the other two dimensions 
are set equal to b, then 

n = 10—^ (cm 4 ) 
ab L 

where M t is given in M , and a and b in parsecs. 

A rough estimate of the age of the shell can be obtained using 
a simple model to describe the expansion of a shell created by 
a continuous injection of mechanical energy or by a supernova 
explosion. In this way, the dynam ical age of GS 10 0-02^41 can 
be estimated as fd yn = a R / V e dWeaver et al.lll977h . where R is 
the radius of the shell (R = yfab) , and a - 0.25 for a radiative 
supernova remnant (SNR) or a = 0.6 for a stellar wind shell. 

Another important parameter that characterizes the shells is 
the kinetic energy, which is given by E k = 0.5 M t V 2 . All the 
relevant parameters of GS 100-02—41 are given in Table [2] 

Even though the CGPS data do not cover the whole area of 
GS 100-02-41, for completeness in Fig. |5] the CGPS data aver- 
aged in the same velocity range than the LAB data are shown. A 
comparison with Fig. [2] clearly shows that the large-scale struc- 
tures are observed in both data sets and, as expected, a lot of 
small scale structures are observed in the high resolution image. 
The dashed ellipse shown in Fig. [5] is the one obtained form the 
fitting of the low resolution data (see Table|2]i. 

3.2. Radio continuum and infrared data. 

Figure [6] shows the observed radio continuum emission at 2695 
MHz (upper panel) and the 60 fim infrared emission (lower 
panel) towards the region where GS 100-02-41 is detected. 

No large scale features that could be interpreted as the coun- 
terpart of the H I structure could be found neither in the radio 
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Table 2. GS 100-02^1 main parameters. 



Parameter 



Value 



Distance (kpc) 


2.8 ± 0.6 


(A>. bo) 


(100! 6, -2 04) 


a 


2!56 


b 


1°69 


e 


-5°7 


a(pc) 


125+ 25 


b(pc) 


83+ 17 


VoCkms- 1 ) 


-41 ± 2 


K e (kms-') 


11 ± 2 


M, (M Q ) 


(1.5± 0.7) xlO 5 


n sh (cm 4 ) 


2.5 ± 0.4 


no (cnT 3 ) 


1.7 + 0.4 


Eu (erg) 


(1.8 ± 0.8) x 10 50 


t dyn (a = 0.25) (Myr) 


2.3 ± 1.1 


t dy „ (a = 0.6) (Myr) 


5.5 ± 2.7 



50 H 7D iO ao 100 




GALACTIC LONGITUDE 



Fig. 5. CGPS averaged HI emission distribution in the veloc- 
ity range from -42.7 to -34.4 km s . The dashed ellipse corre- 
sponds to the ellipse shown in Fig. [2] 



continuum nor in the infrared images. The strong source ob- 
served at both frequencies at (l,b) (102°8, -0P.7) is the well 
known H II re gion S h2-132 , which has been recently studied by 
IVasquezetail fcOlO). Towards the southwest of this H II region, 
a ring nebula related to th e Wolf-Rayet star WR 152 is also ob- 
served (Capp a et al.ll2010t) . On the other hand, three slightly ex- 
tended and much more weaker sources are seen projected onto 
the HI shell. Following the standard nomenclature, they will 
be referred to as G103.39-2.28 (G103, for short), G100.7-0.5 
(G100), and G98.51-1.7 (G98). These three sources will be ana- 
lyzed in Section|4] 

In the infrared image (see Fig. [6] lower panel) two sources, 
IRAS 22036+5306 (IR1) and IRAS 22142+5206 (IR2), located 
at (l,b) = (99?63,-l°85) and (l,b) = (100°38, -3°58), respec- 
tively, are observed. Neither of these sources has a counterpart 
at 2695 MHz. The first source is associated with a post AGB star, 
whilst the second one (IRAS 22142+5206) may be a protostellar 
object dDobashiet al.lll998l) . 

They found IRAS 22142+5206 to be associated with CO 
molecular emission at V = -37.2 kms -1 . 
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Fig. 6. Upper panel. Effelsberg 2695 MHz emission distribu- 
tion in the area of GS 100-02-41. Contour levels are at 65, 250 
and 400 mK. Lower panel. IRIS 60 fin emission distribution. 
Contour levels are at 12, 18 and 24 mJy sr . In both panels the 
ellipse marks the location of the large H I shell. 



3.3. large H I shell origin. 

The origin of GS 100-02—41 may be attributed to the action of 
the stellar wind of several massive stars and their subsequent 
supernova explosions. In the following, based on the derived pa- 
rameters of the shell, we analyze both scenarios. 

An inspection of the Galactic OB Associations in the 
Northe rn Milky Way Galaxy presented by Gar manv & Stencel 
(1992) reveals that several stars belonging to Cep OB 1 are seen 
projected inside the ellipse delineating GS 100-02—41. They are 
ind icated by symbols in Fig. [2] 

iGarmanv & Stencel ( 1 19921) determined for Cep OBI a pho- 
tometric distance modulus of DM = 12.2 mag, which yields 
a distance of 2.75 kpc. On the other hand , based on the stel- 
lar pr oper motions given by Hipparcos, iMel'Nik & Da mbis 
(2009) estimated a distance of 2.8 kpc for this OB association. 
This distance is consistent with the distance estimate derived 
for GS 100-02-41 when non-circular motions are considered 
(lBrand&Biitzlll993l) . To test the possibility that the Cep OBI 
stars lying inside the H I shell could have been capable of creat- 
ing it, an evaluation of the energy that may be injected into the 
ISM by the early type stars of Cep OBI is in place. Is the energy 
provided by these massive stars enough to create GS 100-02- 
41? 

Those stars lying within the boundaries of GS 100-02—41, 
twenty in total, are listed in Table [3] Column 1 gives the 
star's identification, Col. 2 and 3 their galactic coordinates, 
and C ol. 4 their spectral types as given by Garmanv & Stencel 
(fl992h. Based on th e evolutionary track models published by 
Schal ler et al.l (Il992l) . and adoptin g the bolometric magnitude s 
and effective temperatures given by Garmanv & Sten cel (1 19921) . 
we estimated the main sequence (MS) spectral type for each 
star (Col. 5). Column 6 gives an estimate of the star MS life- 
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time (r (MS)) as derived from the stellar models of Schalle retal] Table 4. G98, G100 and G103 parameters. 



( 1992). The values given in this column are a rough estimate as a 
consequence of the uncertainty in the mass adopted for each star. 
Column 7 and 8 give the mass los s rates (M) and wind velocities 
(V w ) taken from lLeifhererl (|T998), respectively. Column 9 gives 
the total wind energy released by each star during its main se- 
quence phase, E w = 0.5 M V^, ?(MS). Given that HDE 235673 
and HDE 235825 are still in the MS, their values are upper lim- 
its. 

Theoretical models predict that only 20 % of th e wind energy 
is con verted into mechanical energy of the shell dWeaver et alJ 
1977 3. In the case of GS 100-02-41 , this implies that for the 
total kinetic energy stored in the shell, E k = (1.8 ± 0.8) x 10 5l) 
erg, a wind energy greater than 9 xlO 50 erg would be required. 
However the analysis of several observed H I shells shows that 
the energy c onversion efficien cy seems to be lower, roughly 
about 2-5 % (ICappa et al.ll2003h . 

From the last column of Table [3] we infer that the total wind 
energy injected during the MS phase of the CepOBl stars is 
~ 29.7 x 10 50 erg, which is enough to create GS 100-02-41 if 
the energy conversion efficiency were ~ 6%. 

It is important to mention that the MS lifetime of all the 
O stars, which are the main energy contributors, is compati- 
ble, within errors, with the dynamical age of the shell (5.5 ± 2.7 
Myr). Moreover, the fact that most of the massive stars located 
inside the shell are evolved stars may explain the absence of ion- 
ized gas related to GS 100-02-41. 

On the other hand, taking into account that the SN rate 
in a typical OB association is about one per 10 5 - 3 x 10 5 
yr. dMcCrav & Kafatosll 19871) . over the lifetime of GS 100-02- 
41 is reasonable to assume that the most massive members of 
Cep OB 1 may have exploded as SN. In this case, the energy of 
the explosions, as well as the energy injected by the stellar winds 
of the SNe progenitors, would have also contributed to the for- 
mation of the shell. There is no detection of a SNR in either 
radio surveys or x-ray surveys, but this could be due to the rel- 
atively large age of GS 100-02—41. However, an evidence that a 
SN explosion could have taken place in the region is the pres- 
ence of the pulsar PSR J2 150+5247, located at (l,b) = (9T.5, 
-0°9) dTavlor et alJl993|). Even though th e estimated distance to 
the pulsar is 5.48 kpc dTavlor et al.ll 19931) . bearing in mind that 
this distance is obtained based on the dispersion measur e and 
a galactic electron density model dTavlor & Cordeslll993h . the 
error is large, and therefore the possibility that the pulsar were 
located at a distance similar to the one of GS 100-02-41 can 
not be ruled out. Were this the case, an upper age limit for the 
progenitor of PSR J2150+5247, assuming it was an early B-type 
star, would be ~3 x 10 7 yr, that is in reasonable agreement with 
the age of the oldest star, HD 209900, belonging to Cep OBI. 



4. Sources G98, G100, and G103 

4.1. Their Nature. 



Source 


G98 


G100 


G103 


Longitude 


98° 5 


100°7 


103° 3 


Latitude 


-V.6 


-a 5 


-2°4 


Size (arcmin) 


73 x 64 


19x23 


54 x 87 


S 4 o8 (Jy) 


2.1 ±0.3 




2.0 ± 0.3 


S1420 (Jy) 


2.0 ± 0.1 


0.21 ± 0.03 


1.50 ±0.15 


S2695 (Jy) 


1.9 ±0.3 


0.19 ±0.03 


1.70 ±0.35 


Spectral index a 


-0.05 ± 0.01 


-0.16 ±0.33 


-0.11 ±0.08 


Sso (Jy) 


410 ± 20 


120 ±6 


1030 ± 50 


S100 (Jy) 


1640 ± 80 


380 ± 20 


1840 ± 90 


Dust temperature (K) 


24.5 ± 0.7 


26.5 ± 0.8 


33.1 ± 1.1 



A clue onto their physical nature (thermal or non-thermal) 
may be provided by the behavior of the flux density as a function 
of frequency . Previous flux d e nsity e stim ations for G1 wer e 
obtained by IWhite & Beckerl d 19921) and iKerton et al] d2007l) . 
IWhite & Becked d!992l) obta i ned a flux density at 1420 MHz of 
456 mJy, while Kerton et~aT1 d2007l) estimated a lower flux den- 
sity, S1420 = 48.8 + 24.0 mJy. To estimate the spectral index 
a (S v oc v - ") of each source, CGPS continuum observations at 
408 and 1420 MHz, and Effelsberg continuum data at 2695 MHz 
were used in this work. 

Flux density determinations, specially for weak sources like 
the ones we are dealing with, are strongly dependent on the 
source angular extent, which in turns depends on the outer- 
most reliable continuum level. Before defining the angular ex- 
tent of each source, the overall galactic continuum emission has 
to be subtracted out. To this end we ap plied the "Back Groun d 
Filtering" method (BGF) developed by ISofue & Reichl d!979l) . 
A l°xl° filtering beam was applied to all sources at each fre- 
quency. 

The derived flux densities are given in Table [4] The quoted 
errors stem from the uncertainty in determining the background 
levels. The spectral index a was derived from a fit to the ob- 
served flux densities, under the assumption that none of these 
sources is optically thick at the lowest frequency (408 MHz). 
For G100 a residual grating ring from the strong nearby Cas A 
source falls on top of the source, making the estimate of the 408 
M Hz flux density mea ningless. The lower flux density obtained 
by iKerton et al.l d2007l) is consistent with the smaller area con- 
sidered for G 1 00 by these authors. On the other hand, the higher 
value given by White & Becker (1992) may be attributed to the 
fact that they used Green Bank data (HPBW ~ 11'). 

Though highly uncertain, the spectral indexes are consistent 
with a thermal nature for G98, G100, and G103; implying that 
they are H II regions. Another argument in favor of this interpre- 
tation comes from the far infrared data. From the measured flux 
de nsities at 60 and 100 /jm a nd following the procedure descibed 
bvlCichowolski et a0 d200lh . the dust temperature was estimated 
for each region (see Table|4]i. The derived temperatures are typ- 
ical of H II regions. The errors quoted for the dust temperatures 
are formal errors from the fitting procedure. 



As mentioned in Section l3T2l three slightly extended sources, la- 
belled G98, G100, and G103, are observed at radio continuum 
and infrared wavelengths, projected onto the border of GS 100- 
02-41. The upper panels of Fig. □ show the CGPS 1420 MHz 
radio continuum image of each source. In these images, discrete 
sources have been removed. The middle panels show the HIRES 
60 /im emission distribution in the area of G98, G100, and G103. 
At both fequencies, to increase the signal-to-noise ratio, the orig- 
inal images were smoothed to a 2-arcmin resolution. 



4.2. Location of G98, G100 and G103. 

We now try to discern whether these ionized regions could be 
physically related to GS 100-02—41. Unfortunately, the radio 
continuum flux densities of G98, G100 and G103 are not strong 
enough to allow us to derive a reasonable H I absorption spec- 
trum in order to attempt to set a limit to their distances. Another 
way to attempt to estimate their distances is to look for signa- 
tures of the interaction among the ionized regions and the HI 



6 



L. A. Suad et al.: GS100-02-41: a new large HI shell in the outer part of the Galaxy. 



Table 3. Cep OBI stars lying inside GS 100-02-^1. 



Star 


/ 


b 


Sp. Type 


MS Sp. Type 


MS lifetime (Myr) 


log (M (M yr-')) 


V»:(kms-') 


£„.(10 M) erg) 


HDE 235673 


98°4 


-r.6 


06.5 V 


06.5 


5.6 


-6.27 


2700 


< 2.2 


HD 209678 


99° 3 


-1!8 


B2I 


BO 


12 


-7.17 


1600 


0.2 


HD 209900 


99! 7 


-1°7 


AO lb 


Bl 


26 


-8.2 


2500 


0.2 


HD 210809 


99° 9 


-3°1 


09 lab 


07 


6.4 


-6.39 


2700 


2.1 


BD +513135 


100°6 


-2°1 


B3II 


B0.5 


14 


-7.3 


1850 


0.2 


HDE 235781 


101° 2 


-?. 6 


B6Ib 


Bl 


22 


-8.2 


2500 


0.2 


BD +53 2820 


101°2 


-1°7 


B0 IV 


BO 


16 


-7.17 


1600 


0.2 


BD +52 2833 


101°4 


-21 


Bl III 


09.5 


10 


-7.04 


2500 


0.5 


BD +53 2837 


101°5 


-2*1 


B2III 


B0.5 


14 


-7.3 


1850 


0.2 


HDE 235783 


101° 7 


-K9 


Bl lb 


B0.5 


14 


-7.3 


1850 


0.2 


BD +53 2843 


101°8 


-T.2 


08III 


07 


6.4 


-6.39 


2700 


2.1 


BD +54 2718 


102°0 


-ff.9 


B2III 


B0.5 


14 


-7.3 


1850 


0.2 


BD +54 2726 


102° 2 


-r.o 


B1.5II 


B0.5 


14 


-7.3 


1850 


0.2 


HDE 2358 13 


102° 4 


-2°.0 


B0 III 


06 


6.3 


-6.14 


2700 


2.6 


HILTNER1106 


102° 5 


-0°7 


B0 III 


08 


6.9 


-6.65 


2600 


0.9 


BD +53 2885 


102° 7 


-T.9 


B2III 


B0.5 


14 


-7.3 


1850 


0.2 


HD 212455 


102° 8 


-V.7 


B6Ib 


B0.5 


14 


-7.3 


1850 


0.2 


HDE 235825 


102° 9 


-1°8 


09 V 


09 


8 


-6.91 


2500 


<0.6 


BD +54 2764 


103°0 


-1°7 


Bl lb 


B0.5 


14 


-7.3 


1850 


0.2 


BD+54 2761 


103° 1 


-1°4 


Q6III 


03* 


4.3 


-5.43 


3200 


<16.3 



References. * : The bolometric magnitude and effective temperature given by Garmany & Stencel ( 1992) lie outside the evolutionary tracks of 
ISchaller et al.l dl992h . We adopted for this star the earliest spectral type. 



and CO gas emission observed around these sources. If the con- 
tinuum sources were related to GS 100-02-41, the HII regions 
would expand within the gas of GS 100-02^-1 and we would 
be able to observe for each source an H I minimum in a velocity 
range compatible with the radial velocity of the large H I shell. 
After a thorough inspection of the CGPS H I data cube, we were 
able to point down HI minima having a good morphological 
correlation with the radio continuum sources. 

Mean brightness temperature of these structures are shown 
in the lower panels of Fig.|7](the original images were smoothed 
to a 2-arcmin resolution). The radial velocity ranges where these 
minima are detected range from -40.2 to -51.8 kms -1 (G98), 
from -50.1 to -58.4 km s 1 (G100), and from -22.1 to -33.6 
kms~' (G103). However, it is important to mention that for the 
averages shown in Fig.|7]only those velocity channels where the 
minima are best defined were considered. To facilitate the com- 
parison among the radio continuum and the H I emission distri- 
butions, the 1420 MHz contour line defining the source's extent 
is superimposed on the H I images. 

The HII regions are clearly surrounded by enhanced HI 
emission. The positional coincidence between the H I features 
and their corresponding H II region suggests that they are phys- 
ically related. It is worth mentioning that, although the entire 
velocity cube was searched for H I features likely to be related 
to the ionized regions, no other peculiar structures but the ones 
shown in Fig. [7] were found. 

We have also analyzed the CO emission distribution in the 
area of G 100 and G103 (unfortunately there is no CO data avail- 
able for G98). We found two molecular clouds probably related 
to G100 in the velocity range form -52.25 to -55.17 kms" 1 
(black contours in Fig. [8}, in coincidence with the velocity range 
where the H I structure was found. Figure[8]also shows the MSX 
Band A (8.28 /mi) image (red colours) and the 5.8 K level at 
1420 MHz (blue contour). Clearly, the 8.28 /mi emission par- 
tially borders the ionized gas. It is worth mentioning that the 
emission observed in the MSX Band A is not detected in the 
other three MSX bands (12.13 /mi, 14.65 //m, and21.3 /mi), sug- 
gesting that the polycyclic aromatic hydrocarbons (PAHs) could 




is - 



Gil*":*!- h'"^^ 



Fig. 8. 8.3 /mi emission distribution in the area of G100. The 
blue contour correspond to the 5.8 K level at 1420 MHz. Black 
contours are associated with the CO emission distribution aver- 
aged between -52.25 and -55.17 kms" 1 , at 0.4, 0.6, 0.8 and 1.0 
K. 



be the main responsible for the emission detected at 8.28 /mi, in- 
dicating the presence of a photo-dissociated region (PDR) in the 
border of G100. An inspection of the CO data cube in the area 
of G 103 do not reveal any CO structure probably related to the 
HII region. 

Bearing in mind that GS 100-02^-1 is detected in the ve- 
locity range from -29 to -51.7 kms -1 and that its baricentral 
velocity is Vo = -41 + 2 kms -1 , we would locate both G98 and 
G100 in the approaching hemisphere of GS 100-02^41, whilst 
G103 would be placed in its receding part. In what fallows we 
adopt for the HII regions the GS 100-02^41 distance. 



4.3. Exciting stars of G98, G100, and G103. 

To identify the exciting stars of G98, G100 , and G103 we have 
inspected the Galactic O Star Catalogue dSota et al.ll2008l) and 
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Fig. 7. Upper panels 1420 MHz radio continuum images in the area of G98 (left), G100 (middle), and G103 (right). Point-like 
sources have been removed. Contour levels are at 6.0, 5.8, and 5.6 K, for G98, G100, and G103, respectively. The black triangle 
symbols indicate the location of the stars HD 235673, ALS 12073, ALS 12071 and ALS 12074 (G98), BD+542684 (G100), and 
ALS 12434, ALS 12443, ALS 12471, ALS 12469 and ALS 12475 (G103). Middle panels: HIRES 60 /xm emission distribution in 
the area of G98, G100 and G103. For comparison with the 1420 MHz emission, contour levels are the same as in upper panels. 
Lower panels: CGPS H I emission distribution averaged in the velocity range from -41.9 to -50.9 km s" 1 in the area of G98 (left), 
from -50.9 to -55.9 kms _I in the area of G100 (middle) and from -29.5 to -31.15 kms 1 in the area of G103 (right). Contour 
levels are the same as for upper and middle panels. In all images the angular resolution is 2-arcmin. 



the OB Star Catalogue (Reed 2003) looking for massive stars 
located in projection toward the H II regions. 

Only one star p rojected within the border of G98 is listed 
bv lSota et alJ (|2008), whilst in the Reed catalogue four stars are 
seen projected towards G98 (on e of them coincides with the one 
catalogued bv lSota et al.1 (120081) 1. one onto G100 , and five onto 
G103. They are listed in Table [5] and their positions are indi- 
cated by triangle symbols in Fig|7] It is worth mentioning that 
HD 235673, the star seen projected onto G98, appears in Table[3] 
as a Cep OB 1 member and therefore was considered as an input 
energy source for the creation of GS 100-02—41. Nonetheless 
its contribution to the overall energy budget is ~ 7 %, so if this 
star were related to G98, its exclusion as being responsible in 
the creation of GS 100-02-41 would not modify the conclusions 
reached in Section [3~3l 

To check whether these stars may have contributed to the 
creation of the observed structures, under the assumption that 
they are located at the same distance than the HII regions (2.8 + 
0.6 kpc), we analyze if their absolute magnitudes are compatible 



with those corresponding to OB-type stars. We need first to es- 
timate the visual absorption in the region. From an averaged H I 
profile (from to -41 km s -1 ) of all over the solid angle covered 
by GS 100-02—41, we derive a total H I column density of about 
Nm = 3 X 10 21 cm" 2 . Using the relation A v = 5.3 x 10~ 22 Ar H i 
(Boh lin et al.ll 19781) . we estimate A v = 1.6 mag. As a check on 
this, for HD 235673 we estimated the visual absorption using 
A v = 3.1 E(B - V), wh ere E(B - V) = (B - V) - (B - V) . 
For ( g-V)n = -0.32 dSchmidt-Kalerl ll982h and (B - V) = 
0.21 ( lHiltner& Johnson! 1 19561) we obtain A „ = 1.6. Both val- 
ues are in agreement w ith the absorption of 1-2 mag given by 
iNeckel & Klard dl98 0). Then, adopting a visual absorption of 1.6 
mag, the absolute magnitudes (see C olumn 5) of the stars listed 
in Table|5]were calc ulated. Based on lMartins & Pled (120061) and 
ISchmidt-Kalerld 19821) . the corresponding spectral types were es- 
timated (see Col. 6 of Table|5j- 

Given that at the assumed distance all the stars listed in Table 
|5]could be O- or early B-type stars, we suggest that they may be 
responsible of creating G98, G100, and G103. 
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Table 5. OB stars probably related to G98, G100, and G103. 



Star 


Galactic coordinates (/, b) 


Sp. Type 


v (mag) 


77T\ 

Ml? 


Sp. Type (o) 


G98 












HD235673 W 


98°36, -L55 


07 


9.14 


-4.7 ± 0.5 


05/08.5 


ALS 12073 


98° 43, -L71 


OB 


12.1 


-1.7 ±0.5 


B2/B5 


ALS 12071 


98°53, -Y.55 


OB 


12.25 


-1.6 ±0.5 


B2/B5 


ALS 12074 


98°55, -1°59 


OB 


12.04 


-1.8 ±0.5 


B2/B5 


G100 












BD +542684 


100°74, -0°55 


OB 


10.8 


-3.0 ± 0.5 


B0/B2 


G103 












ALS 12434 


103° 32, -218 


B5 


10.6 


-3.2 ± 0.5 


B0/B2 


ALS 12443 


103°59, -2°04 


OB 


12.6 


-1.2 ±0.5 


B3/B7 


ALS 12471 


103°59,-2°50 


OB 


12.4 


-1.4 ±0.5 


B2/B7 


ALS 12469 


103°61,-2°44 


B2 


10.6 


-3.3 ± 0.5 


B0/B2 


ALS 12475 


103° 66, -V, 50 


B2 


11.4 


-2.4 ± 0.5 


B1/B3 



References, (a): For the adopted distance of 2.8 ± 0.6 kpc. (b): Star found in both catalogues and member of Cep OB 1 ( Garmanv & Stencel 1992). 



Next, we should consider whether the number of UV ion- 
izing photons needed to keep the continuum sources ionized 
could be provided by their probably exciting stars. The to- 
tal number of Lyman continuum photons is given by N/, vm = 
0.76 x 10 47 T 4 045 v^ z D^ Sy dchaissonlfl976ll . where T 4 is 

the electron temperature in units of 10 4 K, Dk pc is the distance 
in kpc, vghz is the frequency in GHz and S v is the measure 
flux density in Jy. Using the 1420 MHz flux densities given in 
Table|4] a distance of 2.8 kpc and adopting T4 = 1, we obtained 
N Lvm = 1.2 X 10 48 s- 1 for G98, N Lym = 1.3 X 10 47 s" 1 for G100, 
and N Lvm = 1 8 X 10 48 s" 1 for G103. Acc ording to the theoreti- 
cal models of Schaerer & de Koter] (Il997l) . for solar metallicity, 
the estimated N^,,,,, necessary to keep G98, G100, and G103 ion- 
ized, could be provided by a 09.5V, B0. 5 V, and 09.5V stars (or 
two B0V), respectively. This indicates that the join action of the 
stars listed in Table [5] can provide the ionizing photons needed 
for each HII region. In the case of G98, given that the contri- 
bution of HD 235673 is essential to keep this region ionized, we 
conclude that this star is associated with G98 instead of being 
related to the genesis of GS 100-02-41 as it was considered in 
Section [331 This implies that the wind energy provided for this 
star should not have been taken into account for the origin of 
GS 100-02—41. As mentioned before, this fact does not affect 
the scenario proposed for the creation of the large H I shell. 

5. Triggered star formation? 

It is generally believed that expanding sh ells may induce star 
formation at their edges ( Elmegreenl l 1998b . Shells behind shock 
fronts experience gravitational instabilities that may lead to 
the formation of large condensation inside the swept-up gas, 
and some of them may produce new stars. An increasing 
body of observational evidence supports t he importance of 

i ng new sta 

_ _ _ il l2007l ICic 

Cichowolski & Pineauljl201 ll) 

In our case, having an old large H I shell containing several 
HII regions in its edge that seems to be at the same distance, 
we wonder if this could be other case of triggered star forma- 
tion. Were this the case, we should expect an age gradient in 
the region, in the sense that G98, G100, and G103 should be 
younger than GS 100-02—41. In what follows we shall attempt 
to estimate the age of the HII regions G98, G100, and G103, 
and compare them with the dynamical age of GS 100-02^41 es- 
timated in Section l3~T1 (see Table|2|. 



the shell's evolution in creati ng new stars (Patel et al. 1998 
Oev et alJ 120051; lArnal & Cortil 12007; CichowolskietalJ [2009; 



As a rough estimate to the age of eac h region, we evaluat e 
their dynamical ages as td yn = 0.6R/V e dWeaver et al.lll977h . 
The radius R of each region was estimated from the sizes given 
in Table0 We infer/? = 27.8 ± 6.0, 8.5 ± 1.8 and 27.9 + 6.0 
pc, for G98, G100 and G103, respectively, where the errors stem 
from the uncertainty in the distance. Relying on the velocity 
ranges were the H I emission associated with each H II region 
is observed (see Section HJ, expansion velocities of 5.8 kms 1 
(G98), 4.2 km s- 1 (G100) and 5.8 km s _1 (G 103) are assumed. 
Thus, we obtain an age of 2.9 ± 0.9, 1.2 ± 0.5 and 2.9 ± 0.9 Myr 
for G98, G100, and G103, respectively. The quoted errors stem 
from the uncertainties in the sizes and from the assumption that 
the expansion velocities are accurate to within 1.3 kms _I (one 
velocity channel). 

From these estimates, we conclude that G98, G100, and 
G103 are younger than GS 100-02—41. This age difference is 
supported by the fact that most of the stars probably responsible 
of creating GS 100-02-41 have already evolved from the MS, 
while the ionized gas of the H II regions is still observed. 



As mentioned in Section 13.21 in addition to G98, G100, 
and G103, three interesting sources are seen projected onto 
the border of GS 100-02-41: the HII region Sh2-132, the 
ring nebula associated with WR 152, and the infrared source 
IRAS 22142+5206. 

Rega rding Sh2-132 and t he ri ng nebula associate d with 
WR 152, IVasauez et all (MoToh and ICappa et alJ J20Toh found 
neutral gas interacting with the HII regions in the velocity 
ranges from -38 to -53 kms -1 and from —43 to -52 kms _I , 
respectively. For Sh2-132 this is in coincidence with the ion- 
ized velocity g as (Geor gelin & Geqrgelinll 19761 Chu & Tre ffers 



11981b iRevnoldsl 11988b iFich est al.l 1 19901 : lOuireza et al.1 120061. 
Taking into account non-circula r motions in this pa rt of 
the Galaxy dBrand & Blitzl Il993l) . " 
1 et al l d2C~ 



Cappa 1 



Vasau ez et al.1 (2010) and 
(120101) inferred for these regions a kinematical dis- 
tance of 3.5 ± 1.0 kpc. They argued that this value is in close 
agreement with the distance estimates of the main exiting stars 
of the regions, WR 153ab and WR 152. On the other hand, bear- 
ing in mind that GS 100-02—41 is observed in the velocity range 
from -29.0 to -51.7 kms 1 , and that the HII regions are lo- 
cated onto its border, we suggest that both ionized structures may 
be located at the same distance as GS 100-02—41, and that the 
velocity differences are because the HII regions are located in 
the ap proaching part of the large shell. Although Vasauez et al. 
(2010) do not give any age estimate for Sh2-132, the fact that it 
is being ionized by WR 153ab, implies that Sh2-132 should be 
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older than the period of time spent by the progenitor of the WR 
in the main sequence phase. Given that WR 153ab is a WN6.5 
star, the mass of the progenitor was probably a bout 40-50 M© 
dCrowthedl2007h . According to ISchaller et all dl992h the MS 
lifetime for such a star is about 4 Myr. Concerning the ag e of 
the ring nebula associated with WR 152, Cappa et al. (201Q) es- 
timated a dynamical age of 1 Myr for the associated wind blown 
bubble. However, knowing that, as mentioned by ICappa et all 
(2010), large errors are involved in this estimate and bearing in 
mind the time that the progenitor of the WR star spent in the MS 
phase, the age of the H II region is probably higher, closer to the 
one estimated above for Sh2-132. 

With respect to IRAS 22142+5206, as mentioned in Section 
13.21 this infrared source has a massive molecul ar outflow as- 
socia ted with it at the velocity of -37.2 kms 1 (Dobash i et alJ 
1 19981) . The coincidence between this velocity and the veloc- 
ity range where GS 100-02-41 is observed suggests they may 
be located at the same distance. Based on th e total luminos- 
ity of IRAS 22142+5206. lDobashi et ail d 19981) suggest that this 
source will evolve into a late O-type or an early B-type star. 

In summary, the fact that G98, G100, G103, Sh2-132, the 
ring nebula associated with WR 152, and IRAS 22142+5206 lie 
at the edge of GS 100-02-41 and seems to be at the same dis- 
tance, together with the age gradient, suggest that these sources 
could have been triggered by the expansion of GS 100-02—41. 



6. Conclusions 

The large HI shell GS 100-02-41 has been analyzed to study 
the interaction of massive stars with the interstellar medium and, 
in particular, the process of triggered star formation. From this 
analysis we conclude the following: 

1. GS 100-02-41 is a large shell of a radius of about 102 pc 
located at a distance of 2.8 ± 0.6 kpc. The swept up mass 
in the shell is (1.5 ± 0.7) x 10 5 M and the shell density 
n s h = 2.5 ± 0.4 cm" 3 . The shell is expanding at a velocity of 
11 ±2kms~ 1 and its kinetic energy is (1.8+ 0.8) x 10 50 erg. 

2. Several evolved massive stars members of Cep OBI are pro- 
jected inside the large shell. The distance to the OB associa- 
tion is compatible with the kinematical distance of GS 100- 
02^-1 when non-circular motions are considered. An ener- 
getic analysis suggests that the wind energy provided during 
the main sequence phase of the stars could explain the ori- 
gin of the shell. However, taking into account the SN rate in 
OB associations, the energy contribution of a SN explosion 
as well as of its massive progenitor can not be discarded. 

3. From the 2695 MHz radio continuum and 60 /urn infrared im- 
ages, we found three slightly extended sources, labelled G98, 
G100, and G103 projected onto the borders of GS 100-02- 
41 . From the radio flux densities estimated at different wave- 
lengths, the thermal nature of the sources was confirmed by 
the estimation of their spectral indexes. In addition, dust tem- 
peratures were estimated and found to be typical of H II re- 
gions. 

4. An inspection of the 1' CGPS HI data reveals HI minima 
having a good morphological correlation with the HII re- 
gions at velocity ranges compatible with the velocity span- 
ning by GS 100-02-41. This leads to the conclusion that 
G98, G100, and G103 are located at the same distance than 
GS 100-02-41. 

5. From O and OB star catalogues, the massive star candidates 
to be responsible for the ionized gas were identified. 



6. The obtained age difference among the HII regions and the 
shell, together with their relative location leads us to the con- 
clusion that G98, G100, and G103 may have been created as 
a consequence of the action of a strong shock produced by 
the expansion of GS 100-02—41 into the surrounding gas. 
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